Perturbative Approach to the Nonequilibrium Kondo Effect in a Quantum Dot 
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The theory of quantum transport through a dot under a finite bias voltage is developed using 
perturbation theory in the Keldysh formalism. It is found that the Kondo resonance splits into 
double peaks when the voltage exceeds the Kondo temperature, eV > ksTK, which leads to the 
appearance of a second peak in conductance, in addition to the zero-bias peak. The possible relevance 
of the new peak to the 0.7 conductance anomaly observed in quantum point contact is discussed. 
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1. INTRODUCTION 

Ten years after the theoretical predictions of the Kondo 
effect in a quantum dofei^, the zero-bias peak of the 
differential conductance was identified as the Kondo 
resonance^. The observation that the peak value 
reached the unitary limit of 2e 2 /h with decreasing tem- 
perature established unambiguously that the Kondo ef- 
fect is relevant to the transport properties of the quan- 
tum dol£. Clearly the new feature of Kondo transport 
compared with the usual Kondo effect of magnetic impu- 
rities is the nonequilibrium nature, since the current is 
measured with a finite bias voltage. Only in the limit of 
zero bias voltage is the equilibrium condition recovered. 

This theoretical study of the nonequilibrium Kondo ef- 
fect is based on the Keldysh formalism^. In order to treat 
the correlation effect various methods have been used: 
one is perturbation theory with respect to the Coulomb 
interaction U in the dot 7 - and another often used method 
is the non-crossing approximation (NCA) where infinite 
U is assumed^. Concerning the equilibrium Kondo prob- 
lem, it is well known that second order perturbation the- 
ory gives remarkably good results^. For example, the 
density of states with the Kondo resonance in the middle 
of the side peaks corresponding to the energy levels in 
the atomic limit is well described by the theory. How- 
ever, it is not clear whether second order perturbation 
theory works well in nonequilibrium conditions. On the 
other hand, concerning the NCA, it is well known that 
the analyticity is broken in the low temperature limit. 
Thus one cannot discuss the conductance in the unitary 
limit using the NCA. 

One important theoretical issue is that the above men- 
tioned theories predict contradictory results for the den- 
sity of states (DOS) in the nonequilibrium case. The 
NCA predicts splitting of the Kondo resonance under a 
finite bias voltage. On the other hand, in second order 
perturbation theory, the Kondo resonance peak is simply 
suppressed and does not show any particular structure 
in the nonequilibrium situation. Since the Kondo res- 
onance is a manifestation of singlet formation between 
the localized state and the leads, it seems that the dou- 
ble peak structure at chemical potentials of both leads is 



reasonable. Indeed, the double peaks are also obtained 
by other approaches: equations of motion^, real-time di- 
agrammatic formulation 10 - and scaling methods- 11 . More- 
over, it is not clear how the effects of the double peak 
structure, if it is present, would appear in the differential 
conductance. Recently, a proposal was made to measure 
the splitting of the DOS by using three-terminal quantum 
dot o 12 ^ 13 . As a matter of fact, the splitting of the Kondo 
peak was successfully observed by introducing a poten- 
tial difference in the source leadi^. This result is most 
likely due to the splitting of the Kondo resonance, al- 
though the precise geometry is different from the present 
case where a finite voltage is applied between the source 
and drain. It is clearly necessary to study nonequilibrium 
Kondo effects by a better theoretical approach. 

In this paper we analyze the Kondo effect in a quantum 
dot with a finite voltage by using perturbation theory up 
to fourth order in U based on the Keldysh formalism. 
We show that with increasing bias voltage a single Kondo 
peak splits into double peaks at eV ~ UbTk- As a re- 
sult, an anomalous peak of the differential conductance 
appears when eV > IzbTk- At the end of this paper, we 
will also discuss the possible relevance to the experiments 
referred to as the 0.7 conductance anomaly in quantum 
point contacts (QPC)^^^. 



2. MODEL AND CALCULATIONS 

We consider a single-level quantum dot attached to two 
leads. This system is described by the Anderson impurity 
model, 

£ k a cl acr Ckaa- + e d 7l a + U 7l^7l\ 
kaa o 

+ ]T(% QCT cL a d ff +h.c.), (i) 

where a = L, R and CkLa {c-kRcr) annihilates an electron 
in the left(right)-lead, d a annihilates an electron with 
spin a in the dot, and n a = d\d a . The coupling constants 
Vkaa describe the tunneling matrix elements between the 
dot and leads. Nonequilibrium situation is driven by a fi- 



2 



nite difference between and fj, R which are the chemical 
potentials of the leads in both sides. 

For simplicity we concentrate on the symmetric An- 
derson model, where the dot is symmetrically coupled 
with the leads and the energy level of the dot includ- 
ing the Hartree mean field Un a coincides with the cen- 
ter of the potential drop, eV — lil — llr- Then the 
symmetric conditions are stated as Ti, a = r^ CT = T, 
Mi = —l*R = e ^/ 2 and -e d = e d + U. Here T LiRa 
represent resonance width at the chemical potentials, 
^L,Ra{^) = 2vrJ] fc \V k L,Ra\ 2 8(u ~ £fci,fi)- In this paper 
we restrict ourselves to the ground state, T = 0. 

Our first aim is to calculate the density of states in the 
dot, 



estimate the diagonal components, it is more convenient 
to find the Fourier transform of 



1. 



p a {u)) = ImG^(w) 



(2) 



where G r a {u)) is the retarded Green function. Then the 
current through the dot is expressed by 

/OO p p 

duJ r \v - /*H). ( 3 ) 

-OO L(7 I t RcT 

where Jl,r{^>) = 9{^l.r — The differential conduc- 
tance is obtaied from the current by 



G(V) 



dl_ 

dV 



(4) 



3. KELDYSH FORMALISM 
3.1. PERTURBATION THEORY 

Now we briefly sketch perturbation theory based on 
the Keldysh formalism 6 . Each of terms generated by the 
expansion of the S matrix includes integrals along the 
Keldysh contour which starts at t = — oo, passes through 
t = oo and returns back to t = — oo. The branch from 
t = —oo(t = oo) to t = oo(t = — oo) is denoted by the 
index — (+). Thus we need to introduce four types of 
Green functions with these additional indices G a P where 
a and (3 are— or +. 

The interaction terms in the S matrix are given by 
the sum of the Coulomb interaction U at the dot and the 
hybridizations VkL,Ro between the leads and the dot. We 
employ perturbation theory in U based on the Keldysh 
formalism where the Green functions in each order of U 
are renormalized by VkL,Rcr- 

The Green functions of 0th order in U can be explicitly 
evaluated as 

«£-+H = -%2»| 2 Yl r P ,H/ P M, (5) 

p=L,R 

g°+-(u;) = +i\g^)\ 2 £ r p< »(l - / p (o;)),(6) 

p=L,R 

by solving the Dyson equations concerning V kL Ra where 
g°J{u>) = {uj-e d + i/2-{T LrJ {uu)+T Ra {uu))- 1 . In order to 



g°™(t) = 6(-at)g°+-(t) + e(at)g°- + (t), 



(7) 



rather than to solve g® aa directly. 

In the expansion in U there are diagrams in which 
Hartree type self-energies are inserted. Their contribu- 
tions are entirely taken into account by substituting e d —> 
e d +Ufi- a in g® a @ in the diagrams where the Hartree type 
of self-energies are omitted. Using e d + Un^ a = for the 
symmetric case, the new Green functions g^ are 



9a + ( LJ ) = -*2Imfl£(o;)/ e ff(u;), 
g+-(u) = +ahngZ(u)(l-f eB {w)) > 

1 - /effM /eff(w) 



9a M = 

LU 

g+ + (uj) = - 9o 



-iT 



»r' 



(8) 
(9) 

(10) 
(11) 



where Tl jR(T (uj) approximated the values at the chemi- 
cal potentials, /eff(w) = {T La f L (u;) + T Ra f R (uj))/(T LrJ + 
T Ra ) and = (lu + iT)- 1 . 

Then the Dyson equations for U may be written in 
matrix form, 

Gf(u) = g?{u>) + g?{u)lg{u>)G s f{u>). (12) 

Here in each order of £^ 172 it is sufficient to exclude 
Hartree type self-energies. 

Let us evaluate the self-energies by using perturbation 
theory up to the fourth order of U. It is known that 
the perturbative expansion for the Anderson model^ is 
effective and well behaved in the equilibrium case, since 
the exact solution shows a rapid convergence of the per- 
turbation series^. However, in the noncquilibrium situa- 
tion, analysis of perturbation theory has been limited to 
second lowest order. To proceed to higher order calcula- 
tions, the four-point vertex, which is obtained by extend- 
ing the method used by Keldysh for the electron-phonon 
vertes 6 -, is very convenient. This procedure is to insert 
vertices , 

r (0) ai a 2 ,a 3 a4 _ jj( - - _ + + \ 

^ V /Q1Q2 /Q3Q4 /Q1Q2 '03(14/' V 



in each of diagrams, which may be called the Keldysh 
vertices. Here 7";^ = S ai a 2 (7a 2a3 , where (Ta 2 a 3 is tne 
third Pauli matrix. 

All diagrams up to fourth order are shown in Fig. [Tj for 
the symmetric case, where the third-order terms vanish 
in the same way as the equilibrium case. Fig. [IJa) is 
the second-order diagram in which directions of particle 
lines, spins and the Keldysh vertices are depicted. In this 
case, the second order correction of the self-energies is 

£( 2 )7772( w ) = i 2 (_l) f ^^r (0)7172 ' 7374 (14) 

J-00 27T 27T 

gJ^(u Jl )gl 6 J 3 (u J2 )gl 1 J s (uj 1 +u J2 - oj)T^^ . 

We can proceed in the same way to the fourth-order di- 
agrams shown in Fig.[Tjb). Four different representative 
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FIG. 2: The numerically calculated ImEJJ. as a function of u> 
for eV/Y = 0, 0.5 and 2.0. In this example U/T is 6. 



FIG. 1: (a) the second-order diagram and (b) the fourth- 
order diagrams. In (b) all twelve diagrams are generated by 
including the directions of particle lines and spins. Shaded 
squares represent the Keldysh vertices. 



diagrams are illustrated where directions and spin vari- 
able are not specified. By specifying them, one can imme- 
diately obtain all twelve diagrams. Note that Fig. [IJb)- 
1,2,3 and Fig. [IJb)-4 are skeleton and non-skeleton dia- 
grams, respectively. It is tedious but straightforward to 
evaluate the fourth-order contributions, S^ 4 '' 7172 , in the 
same manner as the second order. 

From the self-energy matrix , the imaginary part of Y7 a 
is defined by, 

ImE» = ~ (iS+-H - • (15) 

Then the real part of SJJ. is obtained by the Kramers- 
Kronig relation, 
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ReE;(w) = —V 



Ul' — LU 



did 



(16) 



The Dyson equation for the retarded component is de- 
rived from ea. (fT2")) as 



(17) 



perturbation results correctly reproduce the low-energy 
asymptotic form in the equilibrium limit 9 : 




The nonskeleton diagram does not give any contribu- 
tion to the fourth-order term of the coefficient of (ui/T) 2 . 
Among the skeleton diagrams of the RPA-type diagrams, 
Fig.Htb)-3, give (-3 x 3) for the coefficient of (U/nT) 4 , 
while those from the vertex-correction-type, Fig. []Jb)-l 
and 2, give 3 x (ir 2 — 7), resulting in the small number of 
—3 x (10 — 7r 2 ). The way that cancellation occurs between 
the RPA-type diagrams and the vertex-correction-type 
diagrams changes under a finite bias voltage. 
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3.2. SELF-ENERGIES 

In Fig. [2] the numerically calculated ImS^(w) is shown 
for U/T = 6. In the equilibrium case of eV/T = 0, 
ImE^| e y/ r=0 is dominated by the second-order contribu- 
tion as is well known^, since the fourth-order contribu- 
tions of skeleton diagrams shown in Fig. [ljb)-l,2,3 and 
nonskeleton diagram in Fig. [TJb)-4 almost cancel each 
other except in the low energy limit. In order to an- 
alyze u> <~ 0, it should first be noted that the present 



FIG. 3: The real part of the self-energy for eV/F = 0, 0.5 
and 2.0 with U/T = 6. 

In the vicinity of the equilibrium (for example at 
eV/T = 0.5), -ImS|J(w = 0) becomes slightly larger but 
still keeps the structure of the equilibrium case, which 
means that damping of quasi-particles is simply enhanced 
by real transitions between the leads and the dot. Novel 
features are observed with further increasing the voltage. 
The fourth-order contribution to the coefficient of the 
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(io/T) 2 term changes sign. Even though the second-order 
contribution remains negative, the coefficient of (oj/T) 2 
becomes positive for sufficiently large U/ttT. Therefore, 
the peak of ImEJJ. around w ~ becomes depressed. 
In this case the quasi-particles remain coherent around 
lu = ±eV/2 as is shown for eV/T = 2.0 

Now we turn to ReEJJ. which is shown in Fig. [3J In 
the case of eV/T = 2.0 we see the development of new 
zero-points of Re£^ at u> = ±eV/2 which are absent for 
eV/T = and 0.5. Since the slopes at the zero-points are 
negative, which leads to mass enhancement, the effect of 
damping due to ImEJJ., is reduced. 



4. CONDUCTANCE 

In Fig.[5]the differential conductance defined by eq.(|3]) 
is shown in units of 2e 2 /h for various values of U/T 
as a function of bias voltage. For all U, the zero- 
bias peak starts from the unitary limit. As U is in- 
creased, the width of the zero-bias peak becomes nar- 
rower. For U/T = 4 and 6, a broad peak is seen at 
around eV ~ U. This broad peak corresponds to tun- 
neling processes through the energy levels in the atomic 
limit. Between the zero-bias peak and the broad peak, a 
new peak appears for large U (U/T — 6 in the figure). 



3.3. DENSITY OF STATES 

The density of states is shown in Fig. |4] in units of 
I/71T. In the equilibrium case, the sharp Kondo peak 
develops in the middle of the two peaks at around ±£7/2. 
With a finite voltage, the Kondo peak at u>/T = is 
suppressed but still keeps a single peak with broadening 
for small eV/T, in accordance with the analysis based 
on the Ward identities 19 . When the voltage is further 
increased, the Kondo peak splits into double peaks at u> ~ 
±eV/2, which are located near the chemical potentials of 
the two leads. This behavior is qualitatively consistent 
with the results obtained by the NCA 8 . This change 
occurs when the potential drop eV exceeds fee TV, defined 
as the full-width at half-maximum of the Kondo peak. 
For U/T = 6, IcbTk/T is estimated to be about 0.6. 
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FIG. 4: The density of states in units of I/71T for various 
eV/T with U/T = 6. 

With higher bias the double Kondo peaks located at 
u ~ ±eV/2 merge with the peaks at the atomic limit 
lu ~ ±f//2 and then the latter dominate when eV > U. 
We have checked that this structure coincides with the 
result obtained by the second order perturbation the- 
ory. The origin of this phenomena is that the higher- 
order scatterings become unimportant in the high volt- 
age regime due to the strong dissipative processes be- 
tween the leads and the dot. This result is consistent 
with ref|Tl| in that the nonequilibrium decoherence de- 
stroys the Kondo effect when eV >> ksTx- 
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FIG. 5: The differential conductance in units of 2e 2 /h for 
various values of U/T. The zero-bias peaks reach the unitary 
limit 2e 2 /h. For U/T = 6, the anomalous peak is developed 
in eV > JcbTk in addition to the zero bias peak. 



To understand the origin of the new peak, first we 
rewrite the expression of the conductance, eq.(j4|), as 



G(V) = ^-nT (p a {eV/2) + J 
= G 1 (V)+G 2 (V) 



eV/2 dp^) 



iV/2 



deV 



(18) 



As shown in Fig. O the first term, Gi(V), decreases 
monotonically as V is increased. Concerning the second 
term G ! 2(V r ) first we note that if the range of integra- 
tion integrand is the entire frequency space, the inte- 
gral should vanish due to the sum-rule of the spectral 
weight. For eV < ksTx the spectral weight shifts to 
higher frequencies, giving a negative contribution. When 
eV ~ ksTK a considerable part of positive dp a (uj)/deV 
enters in the integrand, thus giving a less negative con- 
tribution. Therefore we may conclude that the new peak 
appears when the bias voltage exceeds the Kondo energy. 



5. DISCUSSIONS AND SUMMARY 

Finally we discuss the possible relevance of the present 
study to the 0.7 conductance anomaly in QP C 15 ' 16 i 17 . It 
has been suggested that Kondo effect plays significant 
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FIG. 6: The differential conductance for U/F = 6. G(V), 
is given by the sum of Gi(V) and G2(V). The new peak 
originates in the behavior of G2(V). 



role as the origin of the 0.7 structure in QPC because 
the structure is found when the unitary limit of 2e 2 /h^- 
is approached with decreasing temperature. In view of 
the universal nature of the Kondo phenomena, it may be 
interesting to consider experimental results of the QPC 
in light of the present results. Fig. 3 of raffia or Fig. 
1 of refOj] shows that at a gate voltage for which the 
zero-bias conductance reaches the unitary limit the con- 
ductance first drops with increasing the bias voltage and 
then starts to increase again, leading to a second peak at 



finite voltage. The envelope made of these second peaks 
for different gate voltages forms the 0.8 plateau. As the 
temperature is increased the 0.8 plateau is extrapolated 
to the 0.7 structure in the zero-bias limit. In fact, Fig. 
3(b) of reffl7l clearly demonstrates that a second peak 
appears when eV > UbTk- This result seems to support 
the new peak obtained in the present study. 

However the peak-height obtained by the present study 
does not reach to 0.8 as is seen in Fig. [SJ This prob- 
lem may be resolved in the future either by considering 
higher-order corrections or actual level schemes of the 
QPC. 

In summary we have studied the Kondo transport 
through a dot with a finite voltage by using perturba- 
tion theory in the Keldysh formalism. It is shown that 
the splitting of the Kondo resonance occurs when the bias 
voltage exceeds the Kondo temperature. As a result, the 
new peak in the differential conductance appears when 
the electron-electron correlation U is sufficiently strong. 
Finally, we have suggested that the present results are 
relevant also to the 0.7 conductance anomaly. Clearly 
further studies are necessary to elucidate the relation be- 
tween the present analysis and the 0.7 anomaly in the 
QPC. 
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